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ABSTRACT 

Orbital experiments on the behavior of liquid in 
spacecraft are planned. The Sloshsat free-flyer is 
described. Preparation of the experiments, and later 
evaluation, are supported by models of varying 
complexity. The characteristics of the models are 
discussed. Particular attention is given to the 
momentum transfer between the liquid and the 
spacecraft, in connection with the liquid impact that 
may occur at the end of a reorientation maneuver of 
the spacecraft. 

INTRODUCTION 

The European Space Agency (ESA) in conjunction 
with the Netherlands Agency for Aerospace Programs 
has undertaken a project to investigate fluid motion in 
zero-g. ESA has contracted with the Netherlands 
National Aerospace Laboratory (Nationaal Lucht-en 
Ruimtevaartlaboratorium or NLR) for a small free 
flying spacecraft deployed from the shuttle but flying 
in tandem with it. This spacecraft is known as the 
Sloshsat Facility for Liquid Experimentation and 
Verification On-orbit (Sloshsat FLEVO or just 
Sloshsat for short). 1 NASA participates in the 
program and provides the flight opportunity. The 
planned investigations are detailed in an Experiment 
Definition Document. A general objective for all 
experiments is to validate their Computational Fluid 
Dynamics (CFD) models, but NASA has the 
particular goal to demonstrate that propellant can be 
settled at low (Bo < 12) thrust levels, 2 either 
continuously or achieved by on/off modulation. An 
objective of NLR is to validate its Sloshsat Motion 
Simulator (SMS) 3 for use in control algorithms of 
spacecraft. 

nomenclature 

Figure 3 entries: 
m = slug mass 


M = tank mass 
C = tank geometric center 
E = external force point of application 
S = system center of mass(c.o.m.) 

Z = tank c.o.m. 

R = tank radius 

r = re = slug location vector 

v = u + V + flxr = slug inertial velocity 

u = slug relative velocity 

co= slug rotation rate relative to the tank 

V = tank center inertial velocity 
Q = tank rotation rate 

Fl = slug force on tank 
Tl = slug torque on tank 
Fe = thruster force on tank 
Te = thruster torque on tank 

other: 

r = distance of liquid c.o.m. from tank center 
p = density 

V = velocity magnitude 

a = acceleration magnitude 
c = surface tension 
Bo = Bond Number = (2paR~)/o 
We =Weber number = pR V7o 

BACKGROUND 

If a vehicle contains a partially filled tank with liquid, 
problems may arise from the unpredictability of the 
liquid reaction forces in response to a steering 
maneuver. A recent example from space is supplied 
by the NEAR spacecraft that automatically 
interrupted its insertion bum when fuel reaction was 
larger than anticipated. This prevented NEAR from 
orbiting its target asteroid Eros and delayed the 
mission almost a year (The next try was successful in 
putting NEAR in its orbit about asteroid Eros ). 
Another problem, and one that occurs even for 
negligible magnitude of the liquid mass fraction, 
comes from the liquid distribution in a partially filled 
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SYSTEM MODELING 


tank. Most liquid fueled rocket engines require a 
continuous supply of liquid propellant to function. 
Low gravity challenges these engines since the lack of 
separation between liquid and gas will allow gas 
ingestion into an engine, possibly resulting in a 
catastrophic failure of that engine. However, coasting 
in the low gravity environment for long periods of 
time (1/2 hour to days) is required by most upper 
stage rockets to correctly position their payloads. 

Prior to restarting the engines of these rockets, 
something must be done to assure gas free liquid to 
the engines. Satisfactory solutions to liquid slosh and 
fluid management problems would contribute to the 
introduction of high-performance service vehicles 
about ISS. 

DESCRIPTION OF SPACECRAFT 

Sloshsat FLEVO is a free flying spacecraft launched 
from a cross-bay Hitchhiker bridge attached to the 
Space Shuttle. The primary test-bed is an 87-liter 
composite tank filled with 33.5 liters of deionized 
water. Platinum wire rings embedded in the tank wall 
measure capacitance at 270 locations throughout the 
tank. This information is used to determine liquid 
depth over the location to a maximum resolution of 
0.03 m. Miniature thermal dispersion meters provide 
velocity measurements at 10 locations in the tank. A 
cold gas thruster system fed by four gaseous nitrogen 
tanks allows the tank to be maneuvered through a 
broad range of motions. Twelve 0.8 N thrusters 
provide the ability to translate in all three directions 
and rotate about all three axes. Six accelerometers 
and 3 gyroscopes provide detailed information on 
Sloshsat motion. An on-board telemetry system relays 
experiment data back to the Space Shuttle for 
transmission to the ground. Figure 1 shows the 
interior of the Sloshsat. 



Figure 1 Sloshsat Interior 


The investigators with Sloshsat experiments use 
different models according to their objectives. For 
spacecraft dynamics the model must provide reliable 
predictions of the force and torque profile during 
maneuvers. This also implies that the prediction of 
the liquid center of mass (c.o.m. ) should be accurate. 
Since the model is often part of a mission design tool 
it must be simple and CFD codes are generally 
considered not suitable except for final verification 
use. For liquid management investigations the 
emphasis is on good prediction of the liquid shape 
and, again, liquid c.o.m. location. These goals can be 
realized only via CFD. 

The art of model development is disclosed in the 
selection of type and number of variables. The 
number of variables corresponds to the degrees of 
freedom of the model and one generally strives to 
minimize the number. A model needs to be checked 
for correctness, called verification, which is 
accomplished by monitoring conserved quantities 
(momentum, energy) as predicted during operation of 
the model and by recalculation of worked problems 
from the literature. The test of model performance, 
called validation, is realized by comparing its 
predictions with the outcome of physical experiments. 
In some cases comparison with more sophisticated or 
already validated models counts as validation. 

JUSTIFICATION 

Experiments with Sloshsat are deemed necessary for 
the validation of models because the outcome of these 
experiments cannot be predicted with confidence. A 
successful mission is to provide a comprehensive set 
of data to benefit many model developers. 

Previous investigations have indicated that the cause 
of failure to predict may be with two phenomena: 

- the motion of the liquid-solid contact line 

- the damping of the liquid motion. 

SLOSHSAT MODELS 

At present four models are used for Sloshsat mission 
preparation. Two are for spacecraft dynamics analysis 
and use lumped mass liquid, and two are CFD 
programs. The models have been developed 
independently and their predictions are compared for 
verification and for the assessment of range of 
application. In keeping with the Sloshsat mission 
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characteristics the models are meant for large 
amplitude, low Bond number liquid motion. Much of 
the earlier literature covers small amplitude slosh. 
Illustrations of results from the four models are here 
limited to the liquid transfer, or reorientation, 
maneuver. However, many rotational (stability) 
results have been generated also, for preparation of 
Sloshsat control. 

SMS 

The Sloshsat Motion Simulator 3 10 is the first model 
specifically developed for the Sloshsat mission. Its 
architecture is depicted in figure 2; the annotations 
are defined in the Nomenclature. The tank principal 
inertia tensor is also tracked by the model. SMS is 
coded in Fortran and is operated in the EuroS im 
environment. It predicts the 3-D motions of the tank 
and of the liquid c.o.m. The liquid c.o.m. is not 
constrained to a surface but can be anywhere in a 
spherical volume of the cavity of radius R. This 
freedom is a consequence of the liquid model: a 
sphere of constant mass but variable radius, named 
slug'. The slug has variable moment of inertia and 
internal energy from capillary (surface) potential and 
from kinetics. Its resonance frequencies are also 
variable and have values between those of liquid 
drops and bubbles of slug size. 



Figure 2 Architecture of SMS 


CnmFlo 

The ComFlo CFD code has been under development 
since 1995. 9 Comparable to the well-known FLOW- 
3D code, it is based on the Volume of Fluid (VOF) 
method, but it has important advantages when used 
for Sloshsat performance prediction. The code is 
capable of stable simulation of the coupled liquid- 
solid dynamics for any value of the liquid mass 
fraction. Being developed specifically for low Bond 


number flows, it provides accurate tracking of the 
free surface shape and curvature, and without 
problem starts from or terminates to hydrostatic 
equilibrium configurations. The latter feature is 
often mentioned as lacking in FLOW-3 D, due to its 
development as a gravity slosh method. 

Comparison of ComFlo and SMS predictions of 
interaction force and of liquid c.o.m. trajectory 
showed satisfactory agreement . 

ECLIPSE 


a) 5 seconds b) 33.6 seconds 

Two Sloshsat Thrusters Firing Together for a 
0.33 second pulse every second 
(Effective Bo = 4) 

Figure 3 ECLIPSE simulation result 

The ECLIPSE code has a long history of use with 
reorientation problems. 7 ECLIPSE solves the Navier- 
Stokes equation in a radially symmetric formulation 
using a VOF approach. A color function is used to 
track the free surface position. Free surface forces are 
approximated in this version of ECLIPSE by the 
method of Kothe. 8 An illustration of a prediction of a 
transfer maneuver is in figure 3. Heavy black lines 
show the tank wall and free surface location(s). Fine 
lines show vectors of velocity. Since in each image 
the vectors are scaled from the maximum in that 




NASA/TM— 2000-210475 


3 


0.005 0.00720) 


image, they should be viewed for qualitative rather 
than quantitative information. ECLIPSE simulations 
of the Sloshsat settling maneuvers can be found in 
Chato. 2 

Allipx2 

This is a Matlab model of a two-dimensional rigid 
body (the tank) in contact with a movable point mass 
(the liquid) that is constrained to an (approximate) 
elliptical trajectory. Such a model was developed for 
use in the STARDUST mission 1 ' but allips2 has an 
additional interaction. It is force S that acts normal to 
the direction towards the center of the ellipse, to 
model the capillary action for a minimum energy 
configuration in the tank. Such configuration has all 
liquid at one pole of the tank. 

The constitutive equation for S has two parameters, sf 
and sp, as follows: 

S = sf <p' + sp sin2(p /r 

where (p = the angular coordinate of the point mass on 
its trajectory, and <p' its rate of change, or swirl’. 
Parameter sf controls the viscous friction; sp models 
the capillary potential effect. 

Predictions by allips2 have been compared with 
similar predictions generated by the 3-D model, SMS. 
For the comparison the liquid point mass trajectory 
was set to be a circle, and SMS was extended to also 
provide a capillary potential as governed by sp. 

Figure 4 illustrates the different predictions for the 
velocity components of the liquid c.o.m. during a 
reorientation maneuver of Sloshsat. The differences 
can be explained, and comparisons of other dynamic 
variables show similarly good agreement, which 
serves to verify either model. 

Both allips2 and SMS have shown that the liquid 
transfer duration depends strongly on the control gain 
and on the delay time between the tank rate 
measurement and the instance of correction torque 
activation. The ’elasticity’ of the slug in SMS has a 
significant effect on the initial transfer of momentum 
between the slug and the tank. 



Figure 4 Cartesian components of the liquid c.o.m. 
velocity in a spherical/circular cavity 


As a performance test the experimental results of 
Masica have been investigated with the model. For 
liquid in a cylindrical tank, Masica determined at 
what constant acceleration normal to the tank 
centerline the liquid free surface becomes unstable 
and starts to flow. He found that for water (high 
surface tension) the critical Bond number is about 2.5. 
Liquid with low surface tension gave a value 1.12 
and one may suspect some effect of stiction at the 
contact line for the water value. 

The value of the thrust for this run was reduced to 
0.194 N, capable to accelerate the system at 
0.194/(95.3+33.5) = 0.0015 m/s'. Then the Bond 
number, for water, becomes: 

Bo = paR 2 /<y = 1000 * 0.0015 * 0.228 2 / 0.07 = 1.12 

At this, the low value of Masica, the water motion 
should be stable and not break from its initial side of 
the tank, even if stiction effects are not operative. By 
trial and error was found that the stable condition 
requires a magnitude of 0.0057 for sp, value 0.0056 
is unstable. Figure 5 shows the motion of the liquid 
point mass for these values. 

The capillary force parameter sp = 2aR 2 has value 
0.00728 for the Sloshsat data. This means that the 
critical thrust level is predicted to be 0.25 N, or the 
critical Bond number at the reasonable value 1.46. 
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Figure 5. Liquid c.o.m. angular locations, from the 
minor axis of the elliptical trajectory 

The maneuver that generated figure 3 is predicted by 
allips2 as recorded in figure 6. Two curves are shown 
that correspond to slightly different initial conditions 
for the liquid point mass and a consequently different 
transfer trajectory. The dashed line transfer time is 
comparable to the value found by the CFD simulation 
but the solid line transfer takes considerably longer. 
The longer duration transfer (solid line) has the 
masspoint passing far from the tank c.o.m. As in 
Sloshsat, the tank c.o.m. is not at the geometrical 
center of its cavity; it is offset from the major axis of 
the trajectory ellipse in allips2. 
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Figure 6 Transfer predictions by allips2 


Also of interest to the Sloshsat team are the forces 
exerted by the fluid motion and the commanded 
control thruster firings. Allips2 predictions of these 
quantities, for the transfers of figure 6, are shown in 
figures 7 and 8. It is to be noted that the control 
torque activations are reflected in the peaks of the 
liquid-tank interaction forces. 


*Hpe2(200 2 35000 200 0 0.00726) 



Figure 7 Normal forces during the transfers 


4*1062(200 2 35000 200 0 0 00726) 


o.7 r 



#m» In 6*cond> 


Figure 8 Control torque profiles during the transfers 

Allips2 has also been used to estimate the time for 
liquid motion damp during Sloshsat settling 
maneuvers attached to the STS. The tank and thruster 
parameters were given the STS values, and the 
specified operations scenario was implemented. 
When STS primary thrust of 7740 N ( 2 PCRS 
thrusters) is applied along the tank centerline 
direction for 20seconds, it is found that viscous 
damping parameter sf is to have value 0.05 for the 
liquid to stop swirling after 5 minutes. The damping 
of all liquid motion takes considerably longer. If 
thrust is terminated at the right instance, swirl can be 
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absent , or of negligible importance. Figure 9 shows 
the first 30seconds of the liquid c.o.m trajectory 
during a 20seconds primary RCS thruster firing. 
Figure 10 shows the corresponding development of 
the liquid swirl Sf value is 0.05, annotated in the 
figure. Note that the thrust is cut when the swirl has 
near-maximum magnitude; the liquid circles the tank 
in about 4seconds. 



Figure 9 Liquid c.o.m angular location on STS 



VALIDATION QUANTITIES 

The data that are to be produced by Sloshsat are of 
two types: dynamometric (MSS) data and tank data. 
MSS stands for Motion Sensing Subsystem, and the 
data are the output of its six linear accelerometers 
(Allied Signal QA-3000-010) and three gyroscopes 
(LITEF p-FORS). The inertial properties of the 
empty (dry) Sloshsat are measured accurately which 
allows to use this invariable rigid body as a sensitive 
dynamometer for the determination of liquid force 


and torque. Thruster action is calibrated and so can be 
accounted for, but should be absent during precise 
data collection. The tank data are provided by sensors 
inside the tank such as the Coarse Sensor Array of 
platinum wire rings discussed earlier. 

It is asserted that each of the models discussed before 
can be validated from MSS data only. This implies 
that different liquid behaviors will necessarily show 
up in different force and torque time series. Some 
candidate validation quantities will be discussed next. 
If validation cannot be achieved, the tank data will 
provide the main diagnostic tool. Proper evaluation of 
the data will require CFD support since the signal is 
sometimes difficult to interpret. For example, the 
CS A data between zero and saturation are correlated 
not only with local liquid height but also with liquid 
surface inclination and wetted ring fraction. 

The models that calculate the tank (spacecraft) 
dynamics are also able to generate the MSS data and 
so permit direct comparison between prediction and 
actual behavior. Integration of these data allow 
further assessments, e.g. on linear and angular 
trajectories while events like thruster activation or 
liquid transit termination will show their signatures. 
Processing of the data yields the force and torque 
from the liquid interaction, and with an appropriate 
model for the liquid dynamics also the liquid c.o.m. 
trajectory. Evidently, the force and torque must have 
sufficient magnitude and simulations with SMS 
indicate that a force value of 0.2 N might be 
recoverable. 1 1 

The reconstitution of the interaction force is of 
particular interest. Predictions by e.g. SMS or 
ComFlo show typical spikes in the plot vs. time that 
can be removed by averaging. In SMS a spike is 
generated when the slug is at minimum size but still 
has a velocity component directed to the wall. A 
collision then reverts this component at a chosen 
fractional value. Although plausible, such a 
postulated behavior makes the momentum exchange, 
and the loss of kinetic energy (damping) dependent 
on the instantaneous configuration of the liquid and 
tank c.o.m. locations. 14 Engel 15 derives a formula for 
the loss of kinetic energy following impact of a water 
drop on a glass wall. This result can be interpreted to 
give the reconstitution parameter for the collision. 
Alternative assumptions can be made and may yield 
the better choice for agreement with spacecraft data. 
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SUMMARY 

Sloshsat FLEVO provides an opportunity to obtain 
experimental results on liquid behavior during 
maneuvers in the actual low-gravity environment. 

The motional and tank data allow validation of 
various models that have been developed for the 
prediction of these results. 
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